Implications of multiple high-redshift galaxy clusters. 
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To date, 14 high-redshift [z > 1.0) galaxy clusters with mass measurements have been observed, 
spectroscopically confirmed and are reported in the literature. These objects should be exceedingly 
rare in the standard ACDM model. We conservatively approximate the selection functions of these 
clusters' parent surveys, and quantify the tension between the abundances of massive clusters as 
predicted by the standard ACDM model and the observed ones. We alleviate the tension considering 
non-Gaussian primordial perturbations of the local type, characterized by the parameter /nl and 
derive constraints on /nl arising from the mere existence of these clusters. At the 95% confidence 
level, /nl > 467 with cosmological parameters fixed to their most likely WMAP5 values, or /nl ^ 
123 (at 95% confidence) if we marginalize over WMAP5 parameters priors. In combination with /nl 
constraints from Cosmic Microwave Background and halo bias, this determination implies a scale- 
dependence of /nl at ~ 3 o. Given the assumptions made in the analysis, we expect any future 
improvements to the modeling of the non-Gaussian mass function, survey volumes, or selection 
functions to increase the significance of /nl > found here. In order to reconcile these massive, 
high-z clusters with an /nl = 0, their masses would need to be systematically lowered by 1.5 a or 
the as parameter should be ~ 3 o higher than CMB (and large-scale structure) constraints. The 
existence of these objects is a puzzle: it either represents a challenge to the ACDM paradigme or it 
is an indication that the mass estimates of clusters is dramatically more uncertain than we think. 

PACS numbers: cosmology 



INTRODUCTION 

Recent developments in observational hardware and 
observing techniques have enabled the detection of many 
massive, high-redshift clusters [see, e.g. [HI [7J [TT1 155] . 
which seem to create some tension with the abundance 
predictions of the standard ACDM paradigm [551 I5S] . 
Previous work |12[ 1311 135] have examined how the abun- 
dance of high-redshift massive clusters within the ACDM 
model can be enhanced by allowing the primordial fluc- 
tuations, a relic of inflation, to deviate from a Gaussian 
random field. The most basic models of inflation predict 
a scale invariant power spectrum of density pertabations 
<&, described by a Gaussian random field <p. Probes of the 
very early Universe [e.g.. 151] and the Large Scale Struc- 
ture of the late Universe have shown that this description 
is a good approximation to first order. However, any de- 
viations from the slow-roll, single field, adiabatic vacuum 
state inflation (and more complex inflationary models) 
predict deviations from Gaussianity [see e.g., [31 [HI [571 
and refs. therein], which are of interest because they 
1) Modify the number of high-redshift clusters, reliev- 
ing tension between theory and observation, and 2) Al- 
low an observational window into early universe physics. 
The non-Gaussian corrections may be charaterised by the 
coefficient / NL [IE1 HOI EH [57] , which affects the initial 



potential field as 

* = + /nl(0 2 -(0 2 » ■ (1) 

in the so-called local non-Gaussianity case. 

Observations of the Cosmic Microwave Background 
(CMB) WMAP3 by Yadav & Wandelt [BD], measured 
/nl to be within 27 < /nl < 147 (at the 95% confi- 
dence level). More recently, Komatsu & other [28] find 
-10 < /nl < 74 (at 95% C.L.), consistent with the 
above range but also consistent with zero. The CMB 
constrains /nl at large scales (< 0.03/i/Mpc), but on 
smaller scales the Large Scale Structure (LSS) can also 
constrain /nl through the clustering [see e.g. , [49l |56l and 
refs. therein] and abundances of massive halos [see e.g., 
[3TJ[33]. Measurements of /nl using LSS, provides com- 
plementary constraints to the CMB and probes any scale 
dependence of /nl- Considering the scale-dependence on 
halo bias induced by local non-Gaussianity, |59) obtain 
/nl - 53 ± 25 at la, (10 < /nl < 106 at 95% con- 
fidence) from the NVSS survey; this signal comes from 
scales k ~ 0.03/i/Mpc. 

The detection of the high-redshift cluster of galaxies 
XMMUJ2235. 3+2557 g0 and a Hubble Space Telescope 
weak lensing mass measurement [24] , allowed [26] to show 
how the tension between /nl = 0, ACDM (which pre- 
dicts ~ 2 x 10~ 3 such clusters) and this cluster could 



2 



be alleviated with values of 150 < /nl < 260. Massive 
clusters abundance probes /nl on scales corresponding 
to the Lagrangian radius of the halos; k > 0.1 hj Mpc. 

[2"2"] then calculated at which redshift and mass, the 
most massive cluster in the Universe was expected to be 
found, and how this changed with survey volume. They 
also found that XMMUJ2235 . 3+2557 was more than 2a 
away from / NL = 0, ACDM predictions. 

Finally, [10] formally calculated the constraints which 
could be placed on / NL using XMMUJ2235 . 3+2557. They 
computed the probability that the "most massive" clus- 
ter expected within the survey volume had a mass, 1) 
greater than the 68% upper mass estimate of the clus- 
ter, 2) within the 68% upper and lower bounds on the 
mass estimate, and 3) less than the 68% lower bound 
on the clusters mass. They Poisson sampled from these 
abundances to obtain a probability that a cluster with 
the mass of XMMUJ2235 . 3+2557 was the "most massive" 
system. By exploring how values of /nl modified cluster 
abundances [using [35], they placed constraints on /nl to 
be greater than zero at the 95% significance level. We 
note that /nl > is only one possible explanation of the 
existence of high-redshift massive clusters [see, e.g.,[2]. 

The above studies represent the latest results for con- 
straining /nl on (~ lOMpc) cluster scales, and have con- 
centrated on the above single cluster at high-redshift. We 
extend these previous works by exploring the constraints 
on /nl using 14 high-redshift (z > 1.0) spectroscopically 
confirmed galaxy clusters with masses measured in the 
literature. 

The layout of the paper is thus; we begin by re- 
viewing the theoretical form of the cluster mass func- 
tion and the non-Gaussian correction to it, and con- 
tinue by describing the compilation of a high-redshift 
cluster sample. Here we discuss our conservative as- 
sumptions about the selection functions and survey 
volumes. We then describe our analysis and find 
the best fitting cosmological parameters, followed by 
our conclusions and discussions. Throughout the pa- 
per, unless otherwise stated, we assume a flat ACDM 
model with WMAP5 [2T] cosmological parameters (i.e, 
fi m , h, n s , a 8 = 0.28, 0.705, 0.960, 0.812), and quote 
/nl using the LSS convention, e.g. f$l IB ~ /n£ s /1.3 
[see, e.g., US]. 



termines the number of haloes as a function of mass as 
measured within a radius at which the density contrast 
is 180 times the background matter density p m , and has 
the form, 
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where ctm is the rms variation of the density field, 
smoothed on scales M. For a discussion of the minor 
differences between 180/3 m and 200p m mass functions 
see [55] ■ We use the icosmo[5T] package [IS] to calcu- 
late (7m(z), co- moving distances and other cosmology- 
dependent parameters, and use the functional form of / 
[see Equ. B4 of [35] given by 
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Non-Gaussian corrections to the mass function have 
been proposed in the literature [TJl [3TJ G3J [35] , and over 
the mass and redshift ranges considered here, agree to 
within 10% [see Fig. 5&6 of 12 . These corrections 
are typically written as the ratio of the non-Gaussian to 
Gaussian mass functions 1Z, and are, for example, found 
by lineararising the 3-point expansion of the collapse den- 
sity [as in[2I], or by using saddle point approximations to 
non perturbatively account for higher order corrections 
[as in[35] (MVJ), [although, see [33]. We adopt the MVJ 
prescription to describe how the ratio of the non Gaus- 
sian to Gaussian mass functions change as a function of 
/nl 
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where S3 m describes the normalized skewness of the 
smoothed density field, and can be used to define a 



"skewness per /nl unit" as S^m 
given by 
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THE NON-GAUSSIAN CLUSTER MASS 
FUNCTION 

The theoretical cluster mass function was first written 
down by [U] who assumed spherically collapsed halos, 
and was later improved e.g., [50]. Subsequently, large- 
volume, high resolution N-body simulations have been 
performed and mass functions fitting formulae have been 
found [see, e.g., [5J [551 155] ■ We use the spherical over- 
density Gaussian mass function given by [25] . which de- 



where S ec is the critical density for ellipsoidal gravita- 
tional collapse. [58] recently tested these predictions for 
generic non-Gaussianity, using a suite of N-body simula- 
tions, but due to difficulty in computing the initial con- 
ditions, they probed relatively low mass (< 5 x 10 14 M ) 
systems. They found that the MVJ mass function may 
slightly over predict the abundances of massive < 5 x 
10 14 M Q clusters at high-redshift. If this result can be 
extrapolated to more massive clusters at even higher red- 
shifts, then the over prediction of the MVJ non-Gaussian 
mass function will only strengthen the conclusions drawn 
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from this work, as a larger value of /nl will be required 
to fit the observed abundances of massive clusters using 
a more accurate model, implying this analysis is conser- 
vative. 

After the publication of this work, [13] found that the 
exponential fall in the [25] mass function is not enough to 
counter the exponential increase in the non-Gaussian cor- 
rection [35J for very large values of /nl and large masses 
(> 1O 16 M ). They find that the [55] mass function is 
more well behaved for larger values of /nl arid masses, 
but still breaks down at very large scales. We stopped the 
mass function integration at 10 16 M Q just before the [S5J 
mass function breaks down. They additionally checked 
and confirmed the robustness of our method to the choice 
of the mass function and measure a mean value very close 
to that measured here, for the same sample of clusters, 
even after correcting for the mass function approxima- 
tion. In what follows, we only place a lower constraint 
on the value of /nl, and thus our approach is robust to 
the choice of mass function at these lower values of /nl 
and masses considered. 



DATA 



We compile a list of 14 high-redshift (z > 1.0) spectro- 
scopically confirmed clusters with masses measured or 
estimated in the literature, and present them in Table [T] 
We believe this list to represent all known spectroscop- 
ically identified clusters with mass measurements. We 
show the cluster's name, the spectroscopic redshift, the 
cluster mass and mass error converted to M 2 oo (in units 
of 10 14 M Q , assuming an NFW profile [H] if necessary) 
which is the mass enclosed within a radius at which the 
density is 200 times that of the background matter den- 
sity, and the reference to the mass measurement. We 
distinguish clusters detected by X-ray surveys and those 
found using the Sunyaev-Zeldovich [531 hereafter SZ] ef- 
fect. 

Here for each cluster we adopt the mass estimate that 
gave the least tension (best agreement) with /nl = 
ACDM. For an illustrative example consider two cases; 1) 
A cluster mass has a large central value (1 x 1O 15 M0) with 
a large error (4 x 1O 14 M ) , and 2) a cluster has a slightly 
lower mass estimate (7.9 x 1O 14 M ) with a smaller error 
bar 9 x 10 13 M Q ) [see [7]- In our analysis, we find that 
case 1 is more likely to exist in an /nl = ACDM, than 
case 2. Thus, we use case 1 to be conservative. 

We note that mass measurements from different tech- 
niques typicaly agree well, e.g. XMMUJ2235 . 3+2557 had 
mass measurements using weak lensing of 8.3j^'g x 
1O 14 M , and 7.3 ± 1.3 x 1O 14 M gU, and X-ray mass 
measurements of 6 x 10 14 M gTj and 7.7±g \ x 10 14 M 

m 



We also note that potential high-redshift clusters have 
been detected, but not followed up spectroscopically [e.g. 
see [17], so their redshifts, and typically, masses are sub- 
ject to larger uncertainties, if not unknown. This implies 
that our analysis can only place a lower limit on /nl 
as the other clusters may have higher redshifts and/or 
be more massive than the clusters in our sample, which 
would further boost the required value of /nl ■ 

If any of these potential high-redshift clusters candi- 
dates were found to be less massive than those in our 
sample, or at lower redshifts, (and such smaller systems 
are expected in all /nl > 0, ACDM cosmologies), they 
would not detract from these results using the present se- 
lection of clusters, as as our analysis only consider these 
"rare events" that have already been confirmed. 

The ability to detect a cluster, measure its redshift 
and mass for any survey, can be described by the selec- 
tion function. For believable upper and lower limits to be 
placed on cosmological parameters (including /nl) using 
galaxy clusters, the selection function must be under- 
stood. Our analysis uses heterogeneously selected clus- 
ters, so combining the selection functions is non trivial. 
We now describe how we conservatively model the se- 
lection functions for the X-ray and SZ surveys. We note 
that deviations from the conservative modeling, will only 
strengthen our conclusions. 

Selection function 

We split the cluster catalogues into two broad cate- 
gories, those detected using the X-ray by the ROSAT 
and XMM satellites, and those found using the SZ effect 
at South Pole Telescope [9, hereafter SPT]. 

X-ray 

Many of the X-ray surveys have partially overlapping 
footprints, differing flux limits and exposure times. This 
means that some clusters were multiply detected by dis- 
tinct groups, e.g. XMMUJ2235. 3+2557 was originally de- 
tected by the XMM-Newton Distant Cluster Project [40] . 
but was later redetected by the XMM Cluster Survey 
[46] , The combination of all of the X-ray surveys, as per- 
formed here, makes the construction of the full survey 
volume and selection function non-trivial. 

We continue conservatively, by assuming that all X- 
ray surveys had independent footprints (even if they did 
not) and uniform survey volumes (even if some were 
shallower than others), which we choose to be between 
1.0 < z < 2.2 (2.2 represents our estimate of the deepest 
survey limit). We find that our conclusions are stable 
to arbitrary increases of the maximum redshift assumed, 
but will depend on improvements to the modeling of the 
survey footprints and volumes. We reiterate that any 
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TABLE I: We compile a list of high-redshift clusters with mass estimates or measurements from the literature. We show 
the cluster name, redshift, the mass (converted to M200) and la errors, and the mass measurement technique and the mass 
reference. We mark clusters identified from X-ray surveys by + and using the SZ SPT survey by *. The — > indicate the "least 
probable" cluster observed in each of the combined surveys. 



improvements to the conservative selection function and 
footprints adopted here, will make any conclusions drawn 
from this analysis stronger, as a reduced survey volume 
(caused by a smaller footprint or exposure time), or a 
worse selection function (i.e. there are clusters in the 
volume that have not been found) will modify the num- 
ber of observable clusters expected, which will, at best 
not change our results, but at worse, increase tension 
with / NL = ACDM. 

The conservative X-ray survey footprint is 294.5 sq. 
degrees and is composed of 168 sq. degrees from the 
XMM Cluster Survey, 64 sq. deg. from the XMM-Large 
Scale Survey 03], 11 sq. deg. from the XMM-Newton 
Distant Cluster Project, 1.3 sq. deg. from the XMM 
Contiguous survey [2], 17.2 sq. deg. from the Wide 
Angle ROSAT Pointed Survey [32], and 33 sq. deg. from 
the ROSAT Deep survey [15]. 

SZ 

The SZ SPT survey has a well understood selection 
function, and was expected to detect all massive clusters 
above 2 x 10 14 M o [11 US], at all redshifts. We again as- 
sume a survey volume between 1.0 < z < 2.2 and use the 
footprint of 178 sq. degrees. To measure the redshifts 
of clusters detected with the SZ, one needs optical spec- 
troscopic follow up. Not all the identified clusters have 
had their redshifts and masses measured [see [30] j but 
we continue conservatively, by assuming that only clus- 
ters with follow-up were detected. This is conservative 
because future cluster measurements will not relieve the 
tension with /nl = found using the current collection 
of clusters. 




FIG. 1: An Aitoff projection representing the survey foot- 
prints of the combined X-ray survey (shown in green) , the SZ 
SPT survey (yellow), and we also show the SDSS survey foot- 
print (red) for comparison. We note that the X-ray footprints 
are depicted here as a contiguous region, although the actual 
footprints consist of pointings across the whole sky. We rep- 
resent the high-redshift clusters in each survey by crosses and 
triangles. 

Fig [T] is an Aitoff projection representing the survey 
footprints of the combined X-ray survey (shown as a 
green contiguous region, although note that the actual 
X-ray footprint covers many different directions across 
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the full sky), the SZ SPT survey (in yellow), and we also 
show the Sloan Digital Sky Survey [IJ hereafter SDSS] 
survey footprint for comparison (red). The high-redshift 
clusters compiled here, are represented by the crosses and 
triangles. This figure demonstrates how little of the high- 
redshift sky has been observed, and how much volume 
remains to find other potentially massive, high-redshift 
clusters, which may increase the tension with /nl = 
ACDM. 

METHOD AND RESULTS 

Our analysis follows two approaches. First we build 
on the approach of [and refer the reader to §3 of [TO], and 
define the "least probable" (i.e., a combination of most 
massive and highest redshift) clusters in each of the com- 
bined X-ray surveys and the SZ survey, which due to the 
high mass, should also be the easiest to find. We then 
extend the approach of [TU] , by using the existence of the 
compiled cluster sample, including the clusters full mass 
error distributions, to examine the probability that the 
ensemble of clusters could exist in a ACDM universe, and 
probe how the probability increases with /nl- Initially 
we keep the cosmological parameters fixed to WMAP5 
peak values, and then relax this constraint and marginal- 
ize over WMAP5 priors. 

We used the output of both Gaussian (/nl = 0) 
and non-Gaussian (with /nl = 250) N-body simulations 
[obtained from the authors of [55] at a snapshot corre- 
sponding to z = 1.0, to successfully blind test the code 
pipelines. We computed the relative values of /nl needed 
to explain the existence, abundances, and masses of clus- 
ters above > 4x 10 14 M Q , after crudely assuming a survey 
footprint and a redshift slice i.e., a survey geometry. We 
found that at a fixed "probability of existing" , the recov- 
ered value of /nl for the non-Gaussian simulation data 
was always > 225 greater, than that of the Gaussian sim- 
ulation data. For the assumed survey geometry we found 
that the probability of the ensemble of clusters to exist 
was 40% at /nl = in the Gaussian case and < 4% in 
the non-Gaussian case. In the non-Gaussian value 
of /nl = 230 is required to obtain a probability of ex- 
isting to be 40%. We reiterate that the exact recovered 
probability of existence at fixed /nl values, depends on 
the crude conversion of the simulated snapshot volume 
at z = 1.0, to the assumed survey geometry, but the dif- 
ferences between the simulations required a value of /nl 
similar to that inputted into the non-Gaussian simula- 
tions. 

The least probable clusters 

We begin by asking the question, "What is the least 
probable object to be found in each survey assuming 



/nl = 0?". This approach is analogous to determining 
the most massive system in the survey [e.g. [10], but gen- 
eralized to include the redshift-dependence of the mass 
function. 

Assuming the central value for the clusters mass, we 
find that the cluster XMMUJ2235 . 3+2557 is the least prob- 
able X-ray detected object, we expect 5.4 over the full sky 
(and 0.04 in the X-ray survey footprint) at z > 1.39 and 
M > 7.7 x 10 14 M Q using our cosmology and theoretical 
mass function. We also find that SPT-CLJ0546-5345 is 
the least probable SZ detected cluster; we expect only 
12.5 over the full sky (and 0.05 in the survey) with 
M > 1O 15 M and z > 1.06. 

Following [10], we calculate the probability that the 
mass of the "least probable" cluster in each survey falls 
within one of the following three mass bins; 1) less than 
the 1 a mass range of the cluster, 2) within the 1 a mass 
range of the cluster, and 3) greater than the 1 a mass 
range of the cluster. This is accomplished by calculating 
the theoretical cluster abundance within each mass bin, 
and then Poisson sampling from these three abundances 
10 4 times (using the same random number seed for each 
of the three bins), and recording the most massive bin 
which the Poisson samples is > 1. This yields a proba- 
bility that the "most massive" cluster exists is within the 
above mass bins, and within the survey volume. 

We then gradually increase /nl, which boosts the 
abundances of clusters, and Poisson sample from these 
new abundances to re-derive the above probabilities. 
This allows us to place constraints on /nl using the least 
probable observed cluster in each survey. 

In Fig. [2] we show the probability that each ob- 
served massive cluster is the (theoretically predicted) 
"least probable" system in the survey as a function of 
/nl- We note that both clusters provide similar con- 
straints, which, when combined, points to some ten- 
sion with /nl = ACDM. The constraints obtained 
here, are slightly different to that in [TO], due to dif- 
ferences in the assumed survey footprint, mass func- 
tion and cosmological parameters. Note that for ex- 
ample XMMUJ2235. 3+2557 has another (weak lensing- 
based) mass estimate which has a higher central value 
and smaller error-bars. This makes our approach conser- 
vative. 

All clusters 

We proceeded by using the existence, masses and full 
error budgets, of the 14 clusters in the sample. To 
model uncertainties we adopt the following Monte-Carlo 
approach. We log Gaussian random sample from each 
cluster's mass and error 10 4 times producing a set of 
sampled masses Mg, and determine how many clusters 
Ns(M > Ms, z c i us < z < 2.2) one would expect to find 
above each sampled mass and above the redshift of the 
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FIG. 2: The probability that the "least probable" cluster (in 
terms of mass and redshift > z c i us ) in each survey could exist, 
and is the "most massive" cluster in the surveyed volume, as 
a function of /nl- The solid line indicates the probability 
that a cluster more massive than the identified cluster could 
exist, the dotted line is the probability for a cluster within the 
measured 68% mass and error range to be the "most massive 
system", and the dashed line shows the probability that a 
cluster less massive than the cluster is the "most massive" 
system in the survey volume. 



cluster out to edge of survey volume using the mass func- 
tion expression. For each of the 10 4 sampled masses Ms, 
we Poisson sample P° from the predicted abundances 
Nsi and noted if the Poisson sample P°(Ns) > 1, i.e. 
that a cluster more massive than this cluster with a red- 
shift equal to or greater than this cluster could exist. This 
formed a probability Pj, that each cluster i, could ex- 
ist (marginalized over its mass uncertainty), rather than 
forming a probability that the cluster is the "most mas- 
sive", as above i.e., the probability a cluster exists is 
(#P°(iVs) > 1)/10 4 . We then repeated this analysis 
for each of the clusters and multiplied the probabilities 
Pi that each cluster could exist in the surveyed region, 
to produce a combined probability P(/nl) = IIP*, that 
the ensemble of high-redshift clusters could exist in the 
modeled universe. We increased the value of /nl and re- 
peated the analysis to produce a probability distribution 
and stopped the analysis when the P(/nl) = 1? i- e > that 
all the clusters were likely to exist in the cosmological 
model and survey volumes. 

Fig. [3] shows the probability that each cluster could 
exist given the survey volumes and selection function. 
The X-ray and SZ identified clusters are distinguished in 
the figure, but combined in the analysis. We show how 



FIG. 3: The probability that each cluster could exist within 
the survey volume Pi, assuming ACDM with /nl = (black 
symbols). We show X-ray identified clusters by crosses and 
SZ identified clusters by triangles. We also show how the 
probability changes assuming ACDM and setting /nl = 580. 



the probability Pj for each cluster, varies if we change 
/nl from (black symbols) to 580 (red symbols). 

We see that many clusters are unlikely to exist in a 
/nl = 0, ACDM universe, and by multiplying the proba- 
bilities, we find that the probability of the observed Uni- 
verse being well described by this model is 3 x 1CP 3 . 
When /nl = 580 we note that each cluster is more likely 
to exist, and the combined probability = 1 (for our mass 
function), which suggests that this model is a better de- 
scription to the observed Universe [although see, H21 for 
a discussion of the validity of the chosen mass function] . 

In Fig. [4] we plot the combined probability that all the 
clusters could exist as a function of /nl- We see that the 
/nl = model is a poor fit to the observed Universe, 
and by increasing /nl we alleviate tension. We constrain 
467 < /nl at the 95% confidence level using these clus- 
ters. We remind the reader that any improvement in the 
modeling of the survey volumes, footprints or theoretical 
mass function, or the detection of more massive, high- 
redshift clusters, will only increase this result. 



Varying cosmological parameters 

We next simultaneously Gaussian random sample 
from the parameters VLm, ^a, &>k = (1 — Om ~ 
51a), H , er 8 , Wo, n s , ~ 1750 times, using the 
WMAP5 priors (without imposing spatial flatness) and 
record the value of /nl evaluated at P = 0.05, denoted 
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FIG. 4: The probability that the ensemble of clusters could 
exist in a WMAP5 ACDM universe, as a function of /nl ■ 

here as /nl|p(o.05)j which describes the probability of ob- 
serving our 14 clusters P in their surveys P — 0.05 (i.e. 
the exsistence of these clusters in their surveys is allowed 
at 95% C.L.). This procedure is totally analogous to the 
so-called "generalized p-value" for p = 0.05, where the 
uncertainty in the clusters mass and on cosmological pa- 
rameters is effectively marginalized over by treating them 
as "nuisance parameters" with probability distributions 
given by the mass estimates and WMAP constraints. 

In Fig. [5]we show the Id distribution of (generalized) p 
values (so that P > 0.05) as a function of /nl- In other 
words Fig. [5] shows the frequency in our Monte Carlo 
procedure of each value of /nl|p(o.05)- We obtain (123) 
330 < /nl|p(o.05) at 68% (95%) confidence. 

In Fig. [6] we present a selection of two dimensional 
distributions, showing the values of /nl|p(o.05) f° r the 
sampled parameter values, against marginalized distribu- 
tions of; left) the variance of the density field smoothed 
on 8Mpc scales as, and right) the spectral index n s . The 
filled color contours show the 66% (red) and 95% (blue) 
significance levels, and we have marked the peaks in each 
of the distributions by crosses. When viewing these plots, 
one should keep in mind that they represents p value dis- 
tributions for p = 0.05; thus these figures should not 
be interpreted as standard Markov Chain Monte Carlo 
plots. 

We find that /nl is degenerate with a s , but less de- 
generate with all the other varied parameters (we have 
shown only a selection). We can calculate the value of 
as needed for /nl|p(o,05) = by going to lower p values, 
or extrapolating down the line of degeneracy using the 



FIG. 5: The distribution of /nl, which correspond to P(0.05) 
for each Gaussian random sampling of the cosmological 
parameters Om, ^a, tin, fib, Ho, as, Wo, n a using WMAP5 
priors. 

left panel of Fig. [6j resulting in a value of as — 0.87. 
If we only vary as and keep the other parameters fixed 
to their WMAP5 peaks values, we find /nl|p(o.05) = 
when as — 0.89. 

It is interesting to note that Actacama Cosmology 
Telescope found as < 0.86 at 95% CL from upper lim- 
its on the SZ power spectrum |15j . and SPT found 
a s = 0.773 ± 0.025 [32]. The SZ power spectrum sig- 
nal depends very strongly on erg but not as strongly on 
/nl as, for current observations, it is dominated by mas- 
sive (> 1O 14 M ) but lower redshift (z < 1) clusters (see 
[29]). The latest WMAP results alone (combined with 
external data sets) give a more direct, cleaner, measure- 
ment a 8 = 0.801 ± 0.03 (a 3 = 0.809 ± 0.024) [55] ■ The 
high as value necessary to obtain /nl|p(o.05) = is ~ 3a 
away from these constraints. 



CONCLUSIONS AND DISCUSSION 

We compiled a list of 14 high-redshift (z > 1.0) galaxy 
clusters with mass measurements from the literature and 
used their existence to place constraints on the non Gaus- 
sianity parameter /nl- The clusters were identified from 
X-ray surveys and the SZ SPT survey, and we conserva- 
tively assumed a selection function and survey volume. 
We used the theoretical Gaussian mass function of [25 
and the prescription for modifying the cluster abundance 
for non Gaussianities of [35]. We additionally used the 
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FIG. 6: Two dimensional marginalized plots for the value of /nl above which 95% of the probability distribution lies /nl|p(o.05) ; 
against ag and primordial power spectrum spectral index n s . We have represented the peaks in the distributions by thick black 
crosses and the 66% (95%) confidence levels of these p values by red (blue). Note that these figures represent p value distributions 
for p = 0.05, and are not normal MCMC plots. Of all the cosmological parameters explored only as shows a degeneracy with 
/nl. 



output of the Gaussian and non Gaussian N-body simu- 
lations [obtained from the authors of 58 at z = 1.0, to 
successfully blind test the code pipelines. 

We chose to use cluster mass estimates which were per- 
formed assuming a cosmology close to WMAP5 ACDM, 
and to remain conservative, if more than one measure- 
ment technique had been used, we adopted the cluster 
mass and error measurement which allowed for the low- 
est sampled cluster mass. 

We performed two sets of analysis. First we asked the 
question, which is the least probable cluster in each sur- 
vey (this also turns out to be the most massive cluster) 
and asked how likely this cluster was to be the "most mas- 
sive" system in each survey. We found that both massive 
clusters provide some tension with the /nl = WMAP5 
ACDM model, and that by multiplying the probabilities, 
we find that these two clusters have a probability of being 
observed of ~ 30%. 

Using the existence of the 14 clusters, their masses and 
full errors distributions, we then calculated the proba- 
bility that each cluster could exist in the survey. We 
sampled from each cluster's mass and error and calcu- 
lating the expected (Jenkins mass function-predicted) 
abundance above each sampled mass and above the red- 
shift of the cluster, and then Poisson sampled from the 
abundances 62J. We recorded the frequency that the 
Poisson sampled number was greater than or equal to 
one, implying that at least one cluster with the sampled 



mass could exist above the redshift of the cluster in the 
survey volume. We used the frequency of existence to 
construct a probability that each cluster could exist. We 
then combined all probabilities, to obtain a final proba- 
bility that the ensemble of clusters could be found in the 
modeled universe, and we showed how this probability 
changes with /nl- We note that our method allows for 
only a lower limit to be placed on /nl- This is because 
any new clusters, or improvements to the survey volumes, 
or selection functions, will increase tension with /nl = 
ACDM with WMAP priors on cosmological parameters. 

We found that the best fitting models bound /nl 
to be greater than 467 at the 95% confidence level, 
when keeping the WMAP5 parameters fixed at their 
peak values. We also Gaussian random sampled from 
the cosmological parameters CIm, ^a, &k = (1 — — 
17a), ^6, Ho, as, wo, n s using the WMAP5 priors. For 
each realization, we calculated the value of /nl, above 
which 95% of the probability distribution lay. We find 
that the median value of P(0.05) = /nl is 393, and 
drops below /nl = 123, in only ~ 5% of realizations. 
This means that even after marginalizing over cosmolog- 
ical parameters assuming WMAP5 priors, we still find 
/nl|p(o.o5) ^ 123 at the 95% confidence level. 

We have performed several checks: i) the signal is not 
driven by few objects (e.g., only clusters detected in X- 
rays or only those detected in SZ, or only clusters which 
mass estimate is obtained from X-rays etc.) ii) these 
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FIG. 7: A modification of Fig. 8 of Verde [56], with ad- 
ditional /nl measurements from the literature (see text) and 
this work. The colored regions show the scale-dependent mea- 
surements of /nl using the CMB (green colors), the galaxy 
halo bias (blue) and the cluster halo abundances (salmon). 
The enclosed boxes show the 95% confidence levels of each 
measurement. The small x-axis offsets for different measure- 
ments of the same probe is artificial. The lines show the values 
of scale dependent /nl, see text. 



rare events are not evidently clustered in a special patch 
of the sky Hi) cosmological parameters degeneracies: the 
/nl parameter is degenerate only with the erg parame- 
ter. To obtain that /nl = is allowed at 95% C. L., 
the value of eg would have to be ~ 3 a larger than cur- 
rent cosmological (CMB alone and in combination with 
LSS) constraints, iv) all the cluster mass estimates would 
have had to be systematically overestimated by 1.5 a, re- 
gardless of the measurement technique used, to allow the 
ensemble to clusters to be fully compatible with /nl = 
ACDM. 

In Fig. [7J we compare the result obtained here with 
other works, using a modified version of Fig. 8 of [56] . 
We overplot the result on CMB scales (using the at ~ 
0.04/i/Mpc of 27 < /nl < 147, at the 95% confidence 
level by [50] (dark green); of / NL = 32 ± 21 at 1 a by [28] 
(light green); the LSS results at scales ~ 0.4/i/Mpc of 
449 ± 286 at 1 a by [TU] (light salmon, but note that to 
apply an upper constraint, they assume that there will 
be no other clusters found in this footprint as massive 
or more massive than this cluster) ; our result of f^f s > 
123 (so /nl^ 5 > 95, dark salmon); the result using a 
measurement of the non Gaussian scale dependent bias 
at scales - 0.1/i/Mpc of -77 < /nl < 47 at the 95% 
C.L. and peaked at /nl = 8 by [ST] (light blue); and the 



result /nl ~ 53 ± 25 at 1 a (10 < / NL < 106 at the 95% 
C.L.) by [59] (dark blue). They also obtained a similar, 
fully consistent, constraint from the SDSS quasar sample 
(/nl = 58 ± 24). For our application here we use the 
NVSS numbers. 

We used these measurements to constrain the non 
Gaussian spectral index n^Gi defined by [31 , 



/) 



nl 



/nl( 



k* 



(6) 



where the * indicates the CMB pivot scale, k* — 
0.04/i/Mpc. Note that this scale-dependence parame- 
terization does not allow /nl to change sign, so in the 
following approach only /nl > is sampled by our pro- 
cedure. This (theoreticaly-imposed) prior is not too im- 
portant as /nl < for only a small region with relatively 
low probability (recall that [35J finds /nl > at 1.5 a). 

Due to our inability to reliably place an upper con- 
straint on /nl (see the introduction to the data section 
for justification), we assumed a log normal distribution 
for f^ B with a mean of 5.69 and a = 0.212. 

We sampled from the measured values of /nl, while 
keeping k fixed to the central value, and found the 
best fitting curve (using MPFIT[6"3"j) and recorded the 
value of npjG at each pass. The distribution of n^G is 
described by ting = 0.50 ± 0.19 at la, which is a 2.6 a 
detection of scale dependent bias, using [BU], [SJJ and 
our result, or timg = 0.95 ± 0.23 at ler, which is a 4.0 a 
detection of scale dependent bias, using [55], [SS] and our 
result, or ting — 0.93 ± 0.23 at la, using [35], [ST] and 
our result. All of these constraints are in agreement with 
|10) . Since these sets of analysis are not independent, 
the differing results highlight some possible systematics 
effects. We show these lines of best fit on Fig. [7] 



For a non flat distribution of objects, each with an 
observed error, we must account for more objects to be 
scattered into some part of the distribution than are scat- 
tered out. This is described by the Eddington bias, and 
occurs here because the number of expected very massive 
clusters above a mass M, is exponentially smaller than 
the expected number of clusters with mass less than M. 
This could allow lower mass clusters to masquerade as 
higher mass clusters, and potentially cause us to over 
estimate /nl- 

The Eddington bias is estimated to be only a fraction 
of the full la mass error used in this work, and we have 
marginalized over the full mass error distribution and 
have therefore removed any of the Eddington bias effects. 

As a worked example we present the cluster XMMU 
J2235. 3-2557. To calculate the true Eddington bias, 
one should adopt the more robust cluster mass estimate 
not, as we have done here, the more conservative 
one. Typically, the more conservative mass estimate is 
the one with the largest mass error. E.g., {38: states 
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the X-ray mass estimate of XMMU J2235. 3-2557 to 
be 7.7±|;i x 1O 14 M . We find that the statistical 
correction to the mass M, is A In M = 0.48 with 
<7 2 lnAjf = 0.16, and the correction for the Eddington 
bias is AlnM = 0.56, which is indeed higher than 
the la statistical correction (although less than 2a). 
Now, if we instead use the weak lensing mass estimate 
M = 8.5 ± 1.7 x 1O 14 M , of the same cluster, we 
obtain a statistical correction of A In M = 0.2 with 
cr 2 lnM = 0.04, and the corresponding Eddington bias 
correction is AlnM = 0.14. The Eddington bias here 
is therefore 3.5 times smaller than the la statistical er- 
ror of the X-ray estimate, which is that used in this work. 



We conclude with the remarks that we have attempted 
to remain very conservative with our choices of selection 
functions and volumes, with the cluster mass estimates, 
and the modeling of the theoretical non Gaussian cluster 
mass function. Any future improvements in the model- 
ing is expected to strengthen the conclusions of this work; 
if the survey volume decreases, or more clusters are fol- 
lowed up spectroscopically and found to be massive, or 
the theoretical non Gaussian mass function modeling is 
improved, the tension with / NL = WMAP5 ACDM will, 
in all cases, increase. The existence of high-redshift mas- 
sive clusters is a puzzle: it represent a challenge to the 
ACDM paradigme if the clusters mass estimates reported 
in the literature (central values and errors) are taken face 
value. These objects grew too massive too fast com- 
pared to the gravitational instability picture in a ACDM 
paradigm. Alternatively this is an indication that mass 
estimates of high-redshift clusters is dramatically more 
uncertain than currently believed. Weak lensing clusters 
mass estimate is an extremely promising approach to test 
this possibility as [e.g., [34] robust and accurate mass es- 
timates are possible. Such an observational effort would 
help address this "too big, too early" puzzle. 
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